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We study the evolution of a strong field plasma using a quantum Vlasov equation with a non- 
Markovian source term and a simple collision term, and calculate the time dependence of the 
energy- and number-density, and the temperature. The evolution of a plasma produced with RHIC- 
like initial conditions is well described by a low density approximation to the source term. However, 
non-Markovian aspects should be retained to obtain an accurate description of the early stages of 
an LHC-like plasma. 
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At extreme temperature, hadronic matter undergoes a 
phase transition to an equilibrated quark gluon plasma. 
Estimates from phenomenological models and lattice- 
QCD indicate a critical temperature for this transition: 
Tc ^ 170 MeV, which corresponds to an energy density 
of 2-3GeV/fm^. It is hoped that this plasma will be 
produced at the Relativistic Heavy-Ion Collider (RHIC) 
and/or the Large Hadron Collider (LHC) pt]. However, 
little is currently understood about the formation, equi- 
libration and hadronisation of the plasma, and herein we 
focus on dynamical aspects of the pre-equilibrium phase. 

We employ a flux-tube model to explore the creation 
of a strong-field plasma and follow its evolution towards 
equilibrium. The model assumes that the collision of two 
heavy nuclei produces a strong background field and a 
region of high energy density, which decays via pair emis- 
sion. The particles produced in this process are acceler- 
ated by the background field, providing a current and a 
field that opposes the background field. This is the back- 
reaction process, which may result in plasma oscillations. 
For QED in an external field it has been studied via mean 
field methods and using a quantum Vlasov equation with 
a Schwingcr-likc source term [g. Collisions between the 
particles damp the plasma oscillations and are necessary 
to equilibrate the plasma, and their effect has been mod- 
elled in the Vlasov equation approach |^-^ . 

For strong fields, non-Markovian aspects of the parti- 
cle production mechanism are very important M,&-ni|. 
Herein we emphasise this, using a relativistic transport 
equation with the non-Markovian source term derived 
in Refs. [|-0 and explored in Refs. [|l]jl2). We apply 
it with impact energy densities of the scales anticipated 
at RHIC and LHC, and study the time evolution of the 
plasma's properties. 

We model the effect of the nucleus-nucleus collision 
by an external, spatially-homogeneous, time-dependent 
Coulomb-gauge vector potential, which defines the longi- 
tudinal direction: A^ = (0, 0,0, ^(i)). The kinetic equa- 
tion describing fermion production in this external field 
isi 



dfip,t) 
dt 



Sip,t) + C{p,t) 



(1) 



where, in contrast to the Schwingcr production rate, the 
source term here is momentum- and time-dependent: 
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(2) 



E{t) = —dA{t)/dt and e is an electric charge, with the 
dynamical phase and total energy, respectively, 



e(t) = I dt'uj{t') , 

•J — OO 

cu{t) = ^el + (pii - eA{t)y 



(3) 
(4) 



where £± = \/vn?- -V pf^ is the transverse energy, rn^ ~ 
^QCD ~ 0-2 GeV/fm sets a typical scale. 

The second term on the right-hand-side (r.h.s.) of 
Eq. (|]) describes collisions between the particles, and we 
employ a simple and widely-used model 



C{p,t) 



rHp.t) - fip,t) 



(5) 



where r^ is the "relaxation time" and f^"^ is the thermal 
equilibrium distribution functions for fermions: 



rip,t) 



1 



exp[c^(p,t)/T(i)] + l 



(6) 



The temperature profile in Eq. ra) is a priori unknown 
and is determined by requiring that the average energy in 
our ensemble is that of a quasi-equilibrium fermion gas, 
Eqs. (|l5|)-(|l7|). We note that the relaxation time approx- 
imation is only valid under conditions of local "quasi- 
equilibrium" . That is difficult to justify in the presence 
of strong fields. Further, a realistic collision term is likely 
to generate memory effects additional to those present in 
the source term, which are our present focus. The im- 
provement of the collision term is therefore an important 
contemporary challenge. 

The kinetic equation, Eq. p), is non-Markovian for 
two reasons: (i) the source term on the r.h.s. requires 
knowledge of the entire history of the evolution of the 



distribution function from t_oo -^ t; and (ii), even in 
the low density limit (/(i) ~ 0), the integrand is a non- 
local function of time as is apparent in the coherent phase 
oscillation term: cos[Q{t) — Q{t')]. Our kinetic formula- 
tion makes possible a simple and direct connection with 
widely used approximations. 

In the low density limit the source term is independent 
of the distribution function 
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2 uj^{t) 



dt 



,eE{t')ei_ 



u^t 



-^cos[2(e(i)-e(t'))]. (7) 



The kinetic equation, Eq. (M, with the source term of 
Eq. (^ has the general solution 



/dt' exp 
-oo 



(8) 



t' -t 



Slip,t') + 



r±'{p,t') 



The effect of back-reactions on the induced field is ac- 
counted for by solving Maxwell's equation: E{t) = —j{t). 
The total electric field, E{t), is the sum of an external 
field, Eexit), and an internal field, Ein{t). Herein we 
assume that the plasma is initially produced by the ex- 
ternal field, excited by an external current, j^xit), such 
as might represent a heavy ion collision 



Eex{t)^-An[bcos\?{t/b)\ 



(9) 



This model electric field "switches-on" at i ~ — 2& and off 
at i ~ 26, with a maximum magnitude of A^/b at i = 0. 
Once this field has vanished only the induced internal 
field remains to create particles and affect their motion. 
Continued spontaneous production of charged particle 
pairs creates a polarisation current, jpoiit), that depends 
on the particle production rate, S{p,t). Meanwhile the 
motion of the existing particles in the plasma generates 
a conduction current, jcond{t), that depends on their mo- 
mentum distribution, f{p,t). The internal current is the 
sum of these two contributions 



^iny^) — Jin — Jcondy^) Jpoly^) 

The renormalised Maxwell equation is IT 

d^P P\\ - eA{t) 



E,n{t) = -2e 



(27r)3 uj{p,t) 



fip,t) 



(10) 



(11) 



l u;{p,t)df{p,t) 
2uj{p,t) dt 



m 



8w4(p,| _ eA{t)) 



The first two terms on the r.h.s. represent the conduction 
and polarisation current, respectively, while the last term 
arises in regularising the polarisation current. The solu- 
tion of the coupled pair, Eqs. (|^) and (|ll|), yields E{t) 
and f{p, t), and makes evident effects such as plasma os- 
cillations and coUisional damping jl^ . 
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FIG. 1. Time evolution of the energy density for RHIC 
(upper panel) and LHC (lower panel) conditions. 

All bulk thermodynamical properties are expressed in 
terms of f{p, t), and of particular interest herein are: the 
energy density 



eTotAt) = ^E'^{t) + e{t), 



e{t) = 2 



d^p 

(2^ 



L0ip,t)f{p,t). 



(12) 
(13) 



with the pressure p{t) = 1/3 e(i); the particle number 
density 



n{t) 



_fp_ 

(27r)3 



f{p,t): 



(14) 



and the temperature profile: T(i), determined from the 
condition 



n{t) n^iit) ' 
d^p 






(27r)3 

d^p 

(2^ 



oj{p,t)r{p.t), 
rip,t), 



(15) 
(16) 
(17) 



which implements our constraint that the average energy 
in the ensemble is equal to that of a quasi-equilibrium 
fermion gas. 

To solve the coupled equations, we assume a trial form 
for E{t), T{t) and solve for f{p,t) from Eq. (|I|). f(v^) 
so obtained yields an iterated E{t), T{t) from Eqs. ([ll]), 
(p^). The procedure is repeated until seed and iterate 
agree within 0.1%. 
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FIG. 2. Time evolution of the electric field (upper panel) 
and the particle density (low panel) for LHC conditions. 

Our primary goal is to demonstrate that non- 
Markovian effects influence the early stage of plasma evo- 
lution. Therefore we contrast the results obtained in both 
the low-density limit, Eq. (^, and with the full source 
term, Eq. (g). The parameters characterising the exter- 
nal field: b = 0.5/Aqcd « 0.5 fm, A^^^C' ^ 4Aqcd 
and Aq^'^' ~ 20Aqcd, are chosen in order to obtain 
initial field energy densities of the magnitude expected 
in experiments, Fig.|]. We choose a relaxation time 
Tr = I/Aqcd ~ Ifm. 

In Fig. |] we see that for RHIC conditions the energy 
density rises rapidly but, after reaching a maximum, de- 
cays monotonically. In this case the full solution and 
the solution obtained in the low density limit are quan- 
titatively identical. For LHC conditions, with an initial 
energy density twenty-times larger, the situation is dif- 
ferent. The solution obtained in the low density limit 
is only a qualitative guide to the plasma's behaviour and 
plasma oscillations are evident on observable time-scales. 
The deviation between the two curves begins when the 
strength of the external field wanes. 

In Fig. g we present the time evolution of the electric 
field and particle number density for LHC conditions. 
The plasma oscillations are damped on a time scalar char- 
acteristic of the collision time. The evolution obtained in 
the low density limit again provides only a qualitative 
guide to the behaviour of the complete solution. The 
plasma oscillation frequency and amplitude are larger in 
the low density limit because damping is less effective in 
the absence of the Pauli blocking factor: 1— 2/(p, i). This 
same effect is responsible for the overestimate of n{t) in 
the low density limit. 
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FIG. 3. Time evolution of the quasi-equilibrium tempera- 
ture for RHIC and LHC initial conditions. For comparison 
we also plot the Ansatz in Eq. (Ilq). 

Our calculations also yield temperature profiles, which 
are depicted in Fig. ^. The RHIC-like source conditions 
yield an initial temperature r^^^^(i = 0) ^ 0.5 GeV, and 
the temperature decreases monotonically with increasing 
t. The LHC-like source conditions yield an initial tem- 
perature twice as large: T^^'^{t = 0) ■-- 0.9 GeV, and the 
temperature oscillates in tune with the energy density. 
In Ref. |l2| the temperature profile was not determined 
self-consistently, instead an Ansatz was used: 



T{t) - Teg + (T„ 



n. 



'tVtl 



(18) 



with T^g = Aqcd, Tm = 2 Teg, tl = lO/T^, which is also 
plotted in Fig. y. It is evident that the Ansatz provides 
a not unreasonable model of RHIC-like conditions. 

We have solved a quantum Vlasov equation under 
conditions that qualitatively mimic those anticipated at 
RHIC and LHC. Under RHIC conditions the low density 
approximation to the source term provides an accurate 
description of the plasma's early stages. However, the 
non-Markovian features of the source term become im- 
portant when the initial energy density is LHC-like and 
generate effects that are likely to be observable. Follow- 
ing this study we anticipate that the analogue of our non- 
Markovian source term would have a significant effect in 
a non- Abelian transport equation M . 
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